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A  new  family  of  capillary  waves 

By  JEAN-MARC  VANDEN-BROECKf  AND 
JOSEPH  B.  KELLER* 

Courant  Institute  of  Mathematical  Sciences,  New  York  University, 

251  Mercer  Street,  New  York,  NY  10012 

(Received  5  April  1979) 

A  new  family  of  finite-amplitude  periodic  progressive  capillary  waves  is  presented. 

They  occur  on  the  surface  of  a  fluid  of  infinite  depth  in  the  absence  of  gravity.  Each 
pair  of  adjacent  waves  touch  at  one  point  and  enclose  a  bubble  at  pressure  P.  P 
depends  upon  the  wave  steepness  s,  which  is  the  vertical  distance  from  trough  to  crest 
divided  by  the  wavelength.  Previously  Crapper  found  a  family  of  waves  without 
bubbles  for  0  <  s  <  s*  =  0-730.  Our  solutions  occur  for  all  s  >  s**  =  0-663,  with  the 
trough  taken  to  be  the  bottom  of  the  bubble.  As  s-eoo,  the  bubbles  become  long  and 
narrow,  while  the  top  surface  tends  to  a  periodic  array  of  semicircles  in  contact  with 
one  another.  The  solutions  were  obtained  by  formulating  the  problem  as  a  nonlinear 
integral  equation  for  the  free  surface.  By  introducing  a  mesh  and  difference  method, 
we  converted  this  equation  into  a  finite  set  of  nonlinear  algebraic  equations.  These 
equations  were  solved  by  Newton’s  method.  Graphs  and  tables  of  the  results  are 
included.  These  waves  enlarge  the  class  of  phenomena  which  can  occur  in  an  ideal 
fluid,  but  they  do  not  seem  to  have  been  observed. 

1.  Introduction 

Crapper  (1957)  found  a  one-parameter  family  of  explicit  exact  solutions  for  two- 
dimensional  periodic  progressive  capillary  waves  of  finite  amplitude  in  water  of 
infinite  depth.  The  parameter  is  the  steepness  s,  which  is  the  vertical  distance  from 
trough  to  crest  divided  by  the  w-avelength.  His  solutions  are  valid  from  s  =  0,  w-hen  the 
surface  is  flat,  to  s  =  «*  =  0-730,  when  each  two  adjacent  waves  touch  one  another  at 
one  point  (see  figure  1).  Fors  >  «*,  his  formulae  yield  overlapping  waves  with  multiple¬ 
valued  velocities,  so  they  are  not  admissible  as  solutions  of  the  physical  problem.  Thus 
the  wave  of  greatest  steepness  in  this  family  is  that  with  s  =  s*. 

We  have  found  a  new  family  of  waves  of  still  greater  steepness.  They  all  have 
adjacent  waves  touching  at  one  point,  just  like  Crapper’s  highest  wave.  Thus  each 
two  adjacent  waves  enclose  a  region  devoid  of  fluid,  which  we  call  a  bubble.  The 
pressure  P(s)  in  the  bubble  is  an  increasing  function  of  the  steepness,  with  P  =  0  at 
s  =  s*  (see  figure  2).  These  waves  exist  for  all  s  >  s**  =  0-663,  with  P(s)  tending  to  a 
constant  as  s  tends  to  +oo.  Thus,  in  the  interval  s**  <  s  <  «*,  both  these  waves  and 
those  found  by  Crapper  exist,  while,  for  s  >  «*,  the  new  waves  exist  while  Crapper’s  do 
not.  Some  typical  profiles  of  the  waves  and  bubbles  are  shown  in  figure  3. 
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Our  investigation  was  suggested  by  the  work  of  Flaherty,  Keller  &  Rubinow  (1072) 
on  post-buckling  behaviour  of  elastic  tubes  and  rings  with  opposite  sides  in  contact, 
and  by  that  of  Flaherty  &  Keller  (1973)  on  contact  problems  involving  &  buckled 
elastica.  The  method  of  solution  is  numerical,  based  upon  the  formulation  of  the 
problem  as  an  integral  equation  for  the  free  surface,  as  in  Schwartz  &  Vanden-Broeck 
(1979)  and  Vanden-Broeck  &  Schwartz  (1979). 

Kinnersley  (1976)  extended  Crapper’s  exact  result  to  fluids  of  finite  depth.  His 
highest  wave  also  occurred  when  adjacent  waves  touched.  Schwartz  &  Vanden- 
Broeck  (1979)  computed  numerically  capillary-gravity  waves  in  wate  *  of  infinite 
depth.  In  all  the  cases  they  considered  the  solutions  were  limited  by  contact  of  adjacent 
waves.  The  present  method  could  be  used  to  extend  their  solutions,  as  well  as  Kinners¬ 
ley ’s,  to  greater  values  of  the  steepness. 


2.  Formulation  and  solution  without  bubbles 

Let  us  consider  two-dimensional  periodic  capillary  waves  on  the  surface  of  an 
infinitely  deep  fluid.  We  assume  that  they  are  symmetric,  have  wavelength  A,  and 
propagate  with  phase  velocity  c  without  change  of  shape.  Then  we  choose  a  reference 
frame  in  which  the  waves  are  steady,  as  is  the  fluid  motion,  which  is  assumed  to  be  a 
potential  flow.  VVe  also  measure  lengths  in  units  of  A  and  velocities  in  units  of  c,  so  that 
all  variables  become  dimensionless.  Then  we  introduce  rectangular  co-ordinates  ( x,y ) 
with  the  origin  at  a  crest,  with  the  x  axis  parallel  to  the  velocity  at  infinite  depth  and 
with  the  y  axis  directed  out  of  the  fluid. 

The  fluid  flow  can  be  described  by  giving  z  =  x  +  iy  as  an  analytic  function  of  the 
complex  potential /  =  in  the  region  ^  0.  We  choose  the  stream  function  \Jr  to 

vanish  on  the  free  surface,  and  we  set  0  =  0  at  a;  =  y  —  0.  Then  the  free  surface  is  given 
parametrically  by  x(<j>,  0),  y($,  0),  which  we  shall  write  as  x(ip),  y(<j>).  We  denote  the 
complex  velocity  by  q  =  u-  iv,  and  by  our  choice  of  variables  u  tends  to  1  and  v  tends 
to  zero  as  ^  tends  to  —  oo. 

Now  the  dynamical  condition  at  the  free  surface  yields  the  Bernoulli  equation 

qq-i+TK=  0  on  ^  =  0.  (2.1) 


Here  q  is  the  complex  conjugate  of  q,  k(</>)  is  the  curvature  of  the  free  surface  (counted 
positive  when  the  centre  of  curvat  ure  is  on  the  fluid  side  of  the  surface)  and  T  is  the 
dimensionless  parameter 


T  =  2  <r/pAc*. 


(2.2) 


In  (2.2)  p  is  the  fluid  density  and  er  is  the  surface  tension  of  the  surface.  The  choice  of 
the  constant  in  (2. 1 )  implies  that  the  pressure  above  the  surface  equals  that  at  infinite 
depth.  Without  loss  of  generality,  we  can  set  these  pressures  equal  to  zero. 

To  determine  the  free  surface  we  note  that  the  function  +  iy^  —  1  is  analytic  in  the 
half-plane  ^  ^  0  and  vanishes  at  \fr  =  —  oo.  Therefore  its  imaginary  part  on  the  line 
^  =  0  is  the  Hilbert  transform  of  its  real  part: 


d<f>’. 


(2.3) 


The  integral  in  (2.3)  is  to  be  understood  in  the  Cauchy  principal  value  sense. 


A 
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We  now  use  the  assumed  periodicity  and  symmetry  of  the  surface  to  rewrite  (2.3) 
in  the  form 

=1  [x^') -  l][cot 7r(<f>' -<{>)- cot n(<f>' +  <p)]d<f>'.  (2.4) 

V  0 

We  next  rewrite  the  condition  (2.1)  as 

(xl  +  y^-l-Tlx^-y^x^lxl  +  y^-i  =  0.  (2.6) 

Finally  we  note  that  the  steepness  s  of  the  waves  can  be  expressed  in  terms  ofy^(^)  by 
the  formula 


For  a  given  value  of  8,  (2.4)-(2.6)  are  a  set  of  integro-differential  equations  for  the 
functions  x^(<f),  y${$)  and  the  parameter  T .  The  surface  shape  x (<j>),  y(<f>)  can  be  ob¬ 
tained  from  the  solution  by  integration,  and  the  wave  speed  c  can  be  found  from  (2.2). 
The  solution  of  these  equations  should  be  Crapper’s  solution. 

To  solve  these  equations,  we  introduce  the  N  +  1  mesh  points  <j>i  =  (t—  l)/2 N  and 
the  unknowns  x\  —  x^,),  i  =  1,  ...,N  +  1.  Then  we  compute  y\  midway  between 
mesh  points  from  (2.4)  using  the  trapezoidal  rule.  Next  we  compute  y\  and  y\  from  the 
y<+4  by  a  centred  four-point  Lagrange  interpolation  formula  and  a  centred  four-point 
difference  formula,  respectively,  using  the  periodicity  and  symmetry  of  x^.  We  also 
compute  x'  from  the  x'  by  a  four-point  difference  formula.  By  using  these  expressions 
in  (2.5)  at  the  points  <j>it  we  obtain  a  system  of  N  +  1  nonlinear  algebraic  equations. 
Then  by  using  (2.4)  in  (2.6)  we  obtain  another  algebraic  equation.  Thus  we  obtain 
N  +  2  equations  for  the  N  +  2  unknowns  xj,  ...,x's+l,T.  We  solve  these  equations  by 
Newton’s  method.  The  entire  numerical  procedure  follows  closely  that  used  by 
Vanden-Broeck  &  Schwartz  (1979)  and  Schwartz  &  Vanden-Broeck  (1979). 

Some  of  the  results  obtained  by  this  method  for  0  <  s  ^  s*  are  shown  in  table  1  and 
in  figure  1.  They  were  obtained  for  various  values  of  8  with  N  =  40  and  N  =  60.  The 
results  for  these  two  values  of  N,  shown  in  table  1 ,  agree  to  four  decimal  places.  They 
also  agree  well  with  Crapper’s  exact  solution,  which  is  given  in  the  last  column  of  the 
table.  This  is  a  good  check  on  our  numerical  method. 

In  figure  1,  a  few  free  surface  profiles  are  shown  for  different  values  of  sup  to  s  =  0-73, 
at  which  the  surface  touches  itself  and  forms  a  bubble.  For  s  >  0-73  our  numerical 
solution  crosses  itself,  just  as  does  Crapper’s  exact  solution.  Then  the  solution  is  not 
physically  meaningful.  In  Crapper’s  solution,  the  streamlines  of  the  flow  for  a  wave  of 
given  steepness  are  themselves  free  surface  profiles  for  waves  of  smaller  steepness,  so 
they  are  shown  that  way  in  figure  f . 

3.  Waves  with  bubbles 

To  obtain  solutions  for  s  ^  s*  we  require  the  free  surface  to  be  in  contact  with  itself 
at  one  point  within  each  wavelength.  Let  us  denote  by  a  the  <f>  co-ordinate  of  the  contact 
point  in  the  interval  0  <  0  <  J.  Then  by  symmetry  we  have 

x,(a)  =  0,  (3.1) 

x(a)  =  j*  x^(<l>')d<f>'  =  i-  (3.2) 

If  we  adjoin  (3. 1)  and  (3.2)  to  (2.4)-(2.6),  we  have  two  additional  equations  but  only 
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8 

o 

II 

* 

N  =  60 

Crapper 

01 

0-32221 

0-32221 

0-32221 

0-2 

0-33365 

0-33385 

0-33365 

0-3 

0-35188 

0-35188 

0-35188 

0-5 

0-40474 

0-40475 

0-40475 

0-73 

0-48423 

0-48430 

0-48430 

Table  1.  Comparison  of  computed  valuee  of  T  for  N  =  40  and  N  =  60 
with  Crapper’s  exact  value. 


Figure  1.  Computed  free  surface  profiles  for  various  values  of  the  steepness  s.  The  unit  of 
length  is  the  wavelength.  These  results  agree  with  Crapper’s  explicit  solution  for  0  $  «  ^  0-73. 

one  additional  unknown,  a.  Therefore  we  cannot  expect  the  enlarged  set  of  equations  to 
have  a  solution  for  any  value  of  a  other  than  «*,  when  the  solution  is  Crapper’s  highest 
wave,  which  satisfies  (3.1)  and  (3.2)  with  a  =  a*  -  0-17.  This  expectation  is  in  agree¬ 
ment  with  the  results  of  numerical  calculation  which  did  yield  Crapper’s  solution  for 
a  =  «*,  but  which  did  not  yield  any  solution  for  a  >  a*. 

The  physical  reason  why  the  enlarged  problem  does  not  have  a  solution  with  a  >  a* 
is  that  it  requires  the  pressure  P  in  the  bubble  to  be  the  same  as  that  above  the  free 
surface,  i.e.  to  be  P  =  0.  It  is  to  be  expected  that  the  pressure  within  a  bubble  will  have 
some  value  other  than  zero,  which  we  cannot  prescribe.  Therefore  we  shall  modify  the 
free  surface  condition  (2.5)  by  introducing  the  unknown  bubble  pressure  P  into  it  in  the 
interval  a  <  <}>  ^  J,  which  corresponds  to  the  bubble  surface.  Thus  we  replace  (2.5)  by 

(4+yV~'~ 1  *«)(*$ +y$)~*  =  {^p0  *  4  (3  3) 

Here  P  is  measured  in  units  of  pc*. 

Now  we  consider  (2.4),  (2.0),  (3.  J )— (3.3)  as  a  set  of  equations  for  x^),  y^(<j>),  T,  a 
and  P,  for  a  given  value  of  a.  We  have  added  the  two  equations  (3.1)  and  (3.2)  to  the 
problem  of  §  2,  and  added  the  two  unknowns  a  and  P,  so  we  expect  the  new  problem  to 
be  solvable. 
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8 

P 

T 

K 

•S 

0-6635 

-1-6 

0-4727 

3-11 

3  8 

0-7 

-0-28 

0-4801 

3-36 

4-08 

0-73 

0 

0  4843 

3-57 

4-31 

0-75 

013 

0-4863 

3-71 

4-46 

0-85 

0-48 

0-4933 

4-43 

5-24 

105 

0-75 

0-4983 

5-93 

6-83 

1-4 

0-89 

0-4998 

8-61 

9-66 

Table  2.  Computed  values  of  P,  T,  K  and  S  for  various  values  of  «. 


Figure  2.  The  bubble  pressure  P  (in  units  of  pc1)  versus  the  steepness  *  for  the  new  family  of 
waves  with  bubbles.  P  tends  to  unity  as  «  ->•  oo  and  P  =  0  at  a  =  0-73.  The  axis  P  =  0  corres¬ 
ponds  to  Crapper's  solution,  which  is  physically  meaningful  for  IU  i  $  0-73  and  unphysical  for 
a  >  0-73. 


We  shall  solve  this  problem  by  a  numerical  method  similar  to  that  of  the  previous 
section.  We  introduce  the  A,  +  A2  +  1  mesh  points  <pi  defined  by 

<f>i  =  <*(*'-  •)/#!,  *  =  1, ....  A,, 

$i  =  <*+ ($-<*)(*-  Aj  -  1)/A„  i  =  A,+  1 . A,  +  AT,+  1 

We  then  introduce  the  corresponding  unknowns  x\  =  Next  we  discretize 

(2.4)  and  (2.5)  at  the  mesh  points.  In  doing  so,  we  observe  from  (3.3)  that  the  curvature 
of  the  surface  jumps  by  —P/T  at  <p  =  0iVi+,  =  a.  Therefore,  in  interpolating  and 
finite  differencing  at  a  point  on  either  side  of  this  discontinuity,  we  use  only  mesh 
points  on  the  same  side. 

By  using  the  differenced  form  of  (2.4)  as  before  to  eliminate  from  the  differenced 
form  of  (2.5),  we  obtain  a  set  of  A,  + A,+  1  nonlinear  algebraic  equations  for  the 
A,  +  A,+  4  unknowns  z'(,  T,  a  and  P.  We  get  three  more  equations  from  (2.6),  (3.1) 
and  (3.2),  thus  obtaining  A,  +  At  +  4  equations  for  the  same  number  of  unknowns.  This 
system  of  equations  is  solved  by  Newton’s  method. 


(3.4) 
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For  s  slightly  greater  than  s*,  the  iterations  are  started  with  Crapper’s  highest  wave 
as  a  first  approximation.  After  the  iterations  converge,  the  solution  is  used  as  an  initial 
approximation  to  compute  the  solution  for  a  slightly  larger  value  of  a,  and  so  on. 
The  program  was  run  on  a  CDC  6600  computer  with  60  mesh  points. 

4.  Discussion  of  results 

Our  main  numerical  results  are  presented  in  table  2  and  in  figures  2-5.  The  bubble 
pressure  P,  computed  for  various  values  of  the  steepness  s,  was  used  to  draw  the  smooth 
curve  in  figure  2.  This  curve  crosses  the  s  axis  at  s  =  0-730,  where  P  =  0,  and  the  corres¬ 
ponding  solution  is  Crapper’s  highest  wave.  His  solution  for  waves  of  smaller  steepness 
is  represented  by  the  segment  of  the  s  axis  from  s  =  0  to  s  ~  0-730.  The  portion  of  the 
s  axis  with  s  >  0-730  represents  the  non-physical  analytic  continuation  of  his  solution, 
while  the  curve  represents  the  physical  continuation.  Figure  3  gives  the  computed 
free  surface  profiles  for  various  values  of  s,  and  shows  that  the  bubble  size  increases  as 
s  increases. 

We  see  that  the  new  solution  exists  both  for  s  >  0-730,  where  it  is  steeper  than 
Crapper’s  steepest  wave,  and  for  0-663  <  s  <  0-730.  Thus  there  are  two  physically 
admissible  solutions  for  0-663  <  s  <  0-730,  Crapper’s  solution  and  the  new  one.  The 
pressure  P  increases  with  s  along  the  curve  in  figure  2,  so  the  bubble  pressure  exceeds 
the  (atmospheric)  pressure  above  the  fluid  for  s  >  s*  and  is  less  than  that  pressure  for 
s  <  s*.  The  largest  value  of  s  for  which  we  have  computed  the  solution  is  s  =  1-5.  We 
believe  that  we  could  compute  the  solution  for  larger  values  of  s  by  increasing  the 
number  of  mesh  points,  and  that  the  solution  exists  for  all  s  >  0-663. 

As  s  increases,  the  fluid  velocity  q  at  and  near  the  surface  decreases  to  zero  as  ,s-*oo. 
It  then  follows  from  (3.3)  that  the  curvature  of  the  free  surface  tends  to  the  constant 
value  —  T~l  for  0  <  <j>  <  a.  Consequently  this  portion  of  the  free  surface  tends  to  a 
quarter  of  a  circle,  the  diameter  of  which  must  be  equal  to  the  wavelength  1 .  Therefore 
its  curvature  tends  to  -  2,  so  T  tends  to  J.  The  curvature  of  the  bubble,  on  the  other 
hand,  tends  to  zero  at  <j>  =  a,  and  therefore  from  (3.3)  the  pressure  P  tends  to  1. 

These  considerations  indicate  that,  as  s  ->  co,  the  surface  tends  to  a  periodic  array  of 
semicircles  in  contact  with  one  another,  the  bubble  boundaries  tend  to  the  two  sides  of 
the  vertical  lines  beneath  the  contact  points,  and  the  fluid  velocity  at  and  near  the 
surface  tends  to  zero.  This  represents  the  inner  asymptotic  form  of  the  solution  as 
s  ->  oo,  for  y  fixed. 

The  outer  asymptotic  form  is  represented  by  the  flow  beneath  a  periodic  array  of 
bubbles,  each  of  which  extends  upward  to  y  =  +  oo  with  a  width  which  tends  to  zero 
as  y  ->  +  oo.  Matching  of  the  two  asymptotic  solutions  shows  that  the  velocity  q  in  the 
inner  solution  is  proportional  to  e~”’,  which  decreases  rapidly  as  s  increases.  The  bubble 
width  is  also  proportional  to  e*''**'-'1. 

The  bubbles  become  very  long  as  a  increases,  but  adjacent  bubbles  cannot  come 
into  contact  with  one  another.  If  they  did,  they  and  the  free  surface  would  enclose  a 
bounded  region  of  fluid  which  would  therefore  have  to  be  at  rest.  Then  the  adjacent 
bubble  surfaces  would  be  circular  arcs  concave  upward.  They  cannot  be  tangent  to  one 
another  without  being  semicircles,  which  would  require  infinite  curvature  at  the  point 
of  tangency.  This  is  not  allowed  by  the  free  surface  condition. 

As  the  steepness  s  decreases  below  0-730,  P  becomes  more  and  more  negative  and  the 
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Figure  4.  The  dimensionless  parameter  T  =  2<r/pAc*  as  a  function  of  the  steepness  s.  The  lower 
curvo  corresponds  to  Trapper's  solution  and  the  upper  curve  to  the  new  solution.  These  curves 
determine  the  wave  velocity  c  in  terms  of  A  atid  «,  so  they  represent  the  nonlinear  dispersion 
equation  for  the  waves. 

bubble  becomes  smaller.  Figure  3  shows  that  the  curvature  at  0  =  ct_  is  positive.  Thus 
the  bubble  cannot  shrink  to  a  point  as  8  decreases.  If  it  did  the  velocity  at  <j>  =  a  would 
tend  to  infinity.  Then  the  curvature  at  0  =  <x_  would  have  to  be  negative  in  order  to 
satisfy  (3.3). 

We  conclude  that  the  smallest  steepness  8**  for  which  our  family  exists  corresponds 
to  zero  curvature  at  <j>  =  a_.  We  found  numerically  that  «**  =  0-663.  As  we  expected, 
the  numerical  scheme  failed  to  converge  for  s  <  s**. 

In  figure  4  we  present  the  values  of  the  parameter  T  versus  the  steepness  s  for  both 
Crapper’s  solution  and  our  new  family.  The  two  curves  cross  at  s  =  0-730.  We  note  that 
T  tends  to  0-5  as  s  -»■  oo.  This  is  in  agreement  with  the  asymptotic  form  of  the  solution 
as  «->oo,  which  is  described  above. 

The  kinetic  energy  K  and  the  surface  energy  S  can  be  computed  from  our  numerical 
solutions.  Following  Longuetr-Higgins  (1975)  and  Schwartz  &  Vanden-Broeck  (1979), 
we  define  K  and  S  by  the  relations 


K  =  4J!  yd'*>+  Y2  jg  yXd^’ 

(4.1) 

5  =  H,[2j*(^+y'2)*#-l]. 

(4.2) 

A  new  family  of  capillary  waves 
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Figure  5.  The  kiiiotic  energy  K  and  tlio  surface  energy  S,  per  wavelength  (in  units  of  a* /pc*)  as 
functions  of  the  steepness  s.  The  lower  portion  of  each  curve  corresponds  to  Cropper’s  solution, 
and  the  upper  portion  to  tho  now  solution.  The  upper  portions  are  very  nearly  straight  lines. 

The  integrals  in  (4.1 )  and  (4.2)  were  computed  by  the  trapezoidal  rule.  Figure  5  shows 
the  kinetic  energy  K  and  the  surface  energy  8  versus  the  steepness  s  for  Crapper’s 
family  of  solutions  and  for  our  new  family  of  solutions. 
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